Entanglement and majorization in (l+l)-dimensional quantum systems 
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Motivated by the idea of entanglement loss along Renormalization Group flows, analytical ma- 
jorization relations are proven for the ground state of (1 + l)-dimensional conformal field theories. 
For any of these theories, majorization is proven to hold in the spectrum of the reduced density 
matrices in a bipartite system when changing the size L of one of the subsystems. Continuous 
majorization along uniparametric flows is also proven as long as part of the conformal structure 
is preserved under the deformation and some monotonicity conditions hold as well. As particular 
examples of our derivations, we study the cases of the XX, Heisenberg and XY quantum spin chains. 
Our results provide in a rigorous way explicit proves for all the majorization conjectures raised by 
Latorre, Liitken, Rico, Vidal and Kitaev in previous papers on quantum spin chains. 
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I. INTRODUCTION 

In the last few years, the emerging field of quantum 
information science [1] has developed tools and tech- 
niques for the analysis of quantum systems which have 
been proved to be useful in other fields of physics. The 
study of many-body Hamiltonians, quantum phase tran- 
sitions, and the quantum correlations -or entanglement- 
that these systems develop, are examples of this interdis- 
ciplinary research. In fact, understanding entanglement 
has been realized as one of the most challenging and in- 
teresting problems in physics [2], 

Another interesting application of the tools of quan- 
tum information science has been the use of majoriza- 
tion theory [4] in order to analyze the structure present 
in the ground state -also called vacuum- of some models 
along Renormalization Group (RG) flows (for a recent 
review on RG, see [3]). Following this idea, Latorre et 
al. [5] proposed that irreversibility along RG flows may 
be rooted in properties of the vacuum only, without ne- 
cessity of accessing the whole Hamiltonian of the system 
and its excited states. The vacuum of a theory may al- 
ready have enough information in order to envisage irre- 
versibility along RG trajectories. Such an irreversibility 
was casted into the idea of an entanglement loss along 
RG flows, which proceeded in three constructive steps 
for (l+l)-dimcnsional quantum systems: first, due to the 
fact that the central charge of a (l+l)-dimcnsional con- 
formal field theory is in fact a genuine measure of the 
bipartite entanglement present in the ground state of the 
system [6, 7], there is a global loss of entanglement due to 
the c-theorem of Zamolodchikov [8], which assures that 
the value of the central charge at the ultraviolet fixed 
point is bigger or equal than its value at the infrared 
fixed point (cjjv > c ir)\ second, given a splitting of the 
system into two contiguous pieces, there is a monotonic 
loss of entanglement due to the monotonicity numerically 
observed for the entanglement entropy between the two 
subsystems along the flow, decreasing when going away 
from the critical fixed (ultraviolet) point; third, this loss 



of entanglement is seen to be fine-grained, since it follows 
from a very strict set of majorization ordering relations, 
which the eigenvalues of the reduced density matrix of 
the subsystems are numerically seen to perfectly obey. 
This last step motivated the authors of [5] to affirm that 
there was a fine-grained entanglement loss along RG 
flows rooted in properties of the vacuum, at least for 
(l+l)-dimensional quantum systems. In fact, a similar 
fine-grained entanglement loss had already been numeri- 
cally observed by Vidal et al. in [6] , for changes in the size 
of the bipartition described by the corresponding ground- 
state density operators, at conformally-invariant critical 
points. 

In this work, we analytically prove the links between 
conformal field theory (CFT) , RG and entanglement that 
were conjectured in the recent papers [5, 6] for quan- 
tum spin chains. We develop, in the bipartite scenario, 
a detailed and analytical study of the majorization prop- 
erties of the eigenvalue spectrum obtained from the re- 
duced density matrices of the ground state for a variety of 
(1 + l)-dimensional quantum models in the bulk. Our ap- 
proach is based on infinitesimal variations of the param- 
eters defining the model -magnetic fields, anisotropies- 
or deformations in the size of the block L for one of the 
subsystems. We prove in these situations that there are 
strict majorization relations underlying the structure of 
the eigenvalues of the considered reduced density matri- 
ces or, as defined in [5], there is a fine-grained entangle- 
ment loss. The result of our study is presented in terms 
of two theorems. On the one hand, we prove exact con- 
tinuous majorization relations in terms of deformations 
of the size of the block L that is considered. On the 
other hand, we are also able to prove continuous ma- 
jorization relations as a function of the parameters defin- 
ing the model. On top we also provide explicit analytical 
examples for models with a boundary based on previous 
work of Peschcl, Kaulkc and Lcgcza [9]. 

This paper is structured as follows: in sec. II we remem- 
ber the concepts of global, monotonous and fine-grained 
entanglement loss, as defined in [5]. In sec. Ill we ana- 
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lytically prove continuous majorization relations for any 
(l + l)-dimensional CFT when the size of the subsystem L 
is changed, and give an example of a similar situation for 
the case of the XX-model with a boundary. In sec. IV we 
prove continuous majorization relations with respect to 
the flows in parameter space for (l + l)-dimcnsional quan- 
tum systems under perturbations which preserve part of 
the conformal structure of the partition function. Again, 
we support our result with the analysis of a similar sit- 
uation for the Heisenberg and XY quantum spin chains 
with a boundary. Finally, sec.V collects the conclusions 
of our study. We also review in appendix A the definition 
of majorization and provide a lemma which will be used 
in our calculations. 



II. GLOBAL, MONOTONOUS AND 
FINE-GRAINED ENTANGLEMENT LOSS 

Consider the pure ground state (or vacuum) |0) of a 
given physical system which depends on a particular set 
of parameters, and let us perform a bipartition of the 
system into two pieces A and B. The density matrix for 
A, describing all the physical observables accessible to A, 
is given by pa = trs(|Sl) (Q\) -and analogously for B-. In 
this section we will focus our discussion on the density 
matrix for the subsystem A, so we will drop the subindex 
A from our notation. Let us consider a change in one - 
for simplicity- of the parameters on which the resultant 
density matrix depends, say, parameter "t", which can 
be either an original parameter of the system or the size 
of the region A. In other words, we make the change 
p(ti) p(t 2 ), where ti ^ t 2 . In order to simplify even 
more our discussion let us assume that t 2 > t\. We 
wish to understand how this variation of the parameter 
alters the inner structure of the ground state and, in 
particular, how does it modify the entanglement between 
the two parties A and B. Because we are considering 
entanglement at two different points t 2 and t%, we assume 
for simplicity that the entanglement between A and B is 
bigger at the point t\ than at the point t 2 , so we have an 
entanglement loss when going from t\ to t?,. 

Our characterization of this entanglement loss will 
progress through three stages, as in [5], refining at ev- 
ery step the underlying ordering of quantum correla- 
tions. These three stages will be respectively called 
global, monotonous and fine-grained entanglement loss. 

a. Global entanglement loss.- The simplest way to 
quantify the loss of entanglement between A and B when 
going from t\ to t% is by means of the entanglement en- 
tropy S(p(t)) = — ti (p(t) In p(t)). Since at t 2 the two 
parties are less entangled than at ti, we have that 



b. Monotonous entanglement loss.- A more refined 
quantification of entanglement loss can be obtained by 
imposing the monotonicity of the derivative of the en- 
tanglement entropy when varying parameter "i" . That 
is, the condition 



implies a stronger condition in the structure of the 
ground state under deformations of the parameter. This 
monotonic behavior of the entanglement entropy is what 
we shall call monotonous entanglement loss. 

c. Fine-grained entanglement loss.- When 
monotonous entanglement loss holds, we can won- 
der whether, in fact, it is the spectra of the underlying 
reduced density matrix the one that becomes more and 
more ordered as we change the value of the parameter. 
It is then natural to ask if it is possible to characterize 
the reordering of the density matrix eigenvalues along 
the flow beyond the simple entropic inequality discussed 
before and thereby unveil some richer structure. The 
finest notion of reordering when changing the param- 
eter is then given by the monotonic majorization (see 
appendix A) of the eigenvalue distribution along the 
flow. If we call p(t) the vector corresponding to the 
probability distribution of the spectra arising from the 
density operator p(t), then the condition 

p(h)^p(t 2 ), (2.3) 

whenever t<i > t\ will reflect the strongest possible order- 
ing of the ground state along the flow. This is what we 
call fine-grained entanglement loss, and it is fine-grained 
since this condition involves a whole tower of inequalities 
to be simultaneously satisfied (see appendix A) . In what 
follows we will see that this precise majorization condi- 
tion will appear in different circumstances when studying 
(1 + l)-dimensional quantum systems. 



III. FINE-GRAINED ENTANGLEMENT LOSS 
WITH THE SIZE OF THE BLOCK IN 
(l + l)-DIMENSIONAL CFT 

A complete analytical study of majorization relations 
for any (1 + l)-dimensional conformal field theory (with- 
out boundaries 1 ) is presented in the bipartite scenario 
when the size of the considered subsystems changes, i.e., 
under deformations in the interval of the accessible region 



Sipih)) > S(p(t 2 )) , (2.1) 

I assessment between points t 2 and t\. 
This is what we shall call global entanglement loss. theory (BCFT) [17]. 



The case in which boundaries are present in the system must be 
which is a global assessment between points t 2 and t x . considered from the point of view of boundary conformal field 
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for one of the two parties. This size will be represented 
by the length L of the space interval for which we con- 
sider the reduced density matrix pl after tracing out all 
the degrees of freedom corresponding to the rest of the 
universe. Our main result in this section can be casted 
into the following theorem: 



Theorem: pl -i p^ if L > L' for all possible (1 
dimensional CFT. 



1) 



Proof: 

Let Z(t,t) = q- b tr(q iLo+L °' ) ) be the partition func- 
tion of a subsystem of size L on a torus [10], where 
q = e 2lx%T t t = ln with n a positive constant, e 

being an ultraviolet cut-off and b = (c + c)/24 a com- 
bination of the holomorphic and antiholomorphic central 
charges that define the universality class of the model. 
The unnormalized density matrix pl can then be writ- 
ten as pi = q~ b q( Lo+L <>\ since pl can be understood as 
a propagator and (Lq + Lq) is the generator of transla- 
tions in time (dilatations in the conformal plane) [10]. 
Furthermore, we have that 



tr(<j (Lo+Lo) ) = 1 + n iq ai + n 2 q a ' 2 + 



(3.1) 



due to the fact that (Lq + Lq) is diagonal in terms of 
highest- weight states \h, h): (Lq + L )\h, h) = (h + 
h)\h, h), with h > and h > 0; the coefficients 
ai, ct2, ■ ■ ■ > 0, Oj+i > cti Vi are related with the scaling 
dimensions of the descendant operators, and ni,U2,--- 
are degeneracies. The normalized distinct eigenvalues of 
Pl are then given by 



Ai 



A 2 = 



1 


(1- 


h niq ai + n 2 q a2 + 


•••) 








(1- 


hni9 Q i +n 2 q a2 + 


...) 




q a d-i) 




(1- 


h nxq^ + n 2 q a2 + 


■■•) 



(3.2) 



Let us define Z(q) = tr(g( Ln+L °)) = (1 + nxq 011 + 
n 2 q a2 + •••). The behavior of the eigenvalues in terms of 
deformations with respect to parameter L follows from, 



dZ(q) _ Z(q) - 1 dq d\n(L/e) 



dL 

and therefore 



dln(L/e) dL 



>0, 



dAi 
dL 



d 

dL 



Z(q) 



< 0. 



(3.3) 



(3.4) 



Because Ai is always the biggest eigenvalue VX, the first 
cumulant automatically satisfies continuous majorization 
when decreasing the size of the interval L. The variation 



of the rest of the eigenvalues A/, I > 1, with respect to L 
reads as follows: 




Z(q) 



Z(q) 




(3.5) 



Let us focus on the second eigenvalue X 2 . Clearly two 
different situations can happen: 



• if (ai - Zi z 9) iq) 1 ) > 0, then since a (i _i) > ai VZ > 
2, we have that (a^ - ^§^f) > V/ > 2, 



which in turn implies that > V/ > 2. From 
this we have that the second cumulant satisfies 



d(Ai + A 2 ) 
dL 



< 



(3.6) 



d>2 



thus fulfilling majorization. The same conclusion 
extends easily in this case to all the remaining cu- 
mulants, and therefore majorization is satisfied by 
the whole probability distribution. 
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d(A!+A 2 ) 



Z(q)-1 
Z(q) 



< 0. then 



dA 2 
dL 



< 0, and therefore 



— — < 0, so the second cumulant satisfies ma- 
jorization, but nothing can be said from this about 
the rest of the remaining cumulants. 

Proceeding with this analysis for each one of the eigen- 
values we see that, if these are monotonically decreasing 
functions of L then majorization is fulfilled for the par- 
ticular cumulant under study, but since o^+i > at Vz 
we notice that once the first monotonically increasing 
eigenvalue is found, majorization is directly satisfied by 
the whole distribution of eigenvalues, therefore pl -< pv 
if L > L' . This proof is valid for all possible (1 + 1)- 
dimensional conformal field theories since it is based only 
on completely general assumptions. □ 



A. Analytical finite-L majorization for the critical 
quantum X-model with a boundary 

Let us give an example of a similar situation to the 
one presented in the previous theorem for the case of 
the quantum XX-model with a boundary, for which the 
exact spectrum of pl can be explicitly computed. The 
Hamiltonian of the model without magnetic field, is given 
by the expression 



(3.7) 



The system described by the XX-model is critical since 
it has no mass gap. Taking the ground state and tracing 
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out all but a block of L contiguous spins, the density 
matrix p^ describing this block can be written, in the 
large L limit, as a thermal state of free fermions (see [9]): 



IV. FINE-GRAINED ENTANGLEMENT LOSS 
ALONG UNIPARAMETRIC FLOWS IN 
(1 + l)-DIMENSIONAL QUANTUM SYSTEMS 



-H' 



PL 



z, 



(3.8) 



Zl being the partition function for a given L, H' = 
X)fc=o e kd\dk, with fermionic operators d\, dk and dis- 
persion relation 



We study in this section strict continuous majorization 
relations along uniparamctric flows, under the conditions 
of integrable deformations and monotonicity of the eigen- 
values in parameter space. The main result of this section 
can be casted into the next theorem: 



2 InL 



(2k + 1) k = 0, 1,. ..,L- 1 . 



(3.9) 



The eigenvalues of the density matrix p^ can then be 
written in terms of non-interactive fermionic modes 



PL(n ,ni, . . . ,n L _i) = — e 



r-v L — l 



pL(n ) ■ • -Pi(rii-i) 



(3.10) 



is the 



with p(n a ) = ie-"° £ °, where Z% = (1 + e 
partition function for the mode a, and n a = 0, 1, Va. It 
is worth noticing that the partition function of the whole 
block Zl can then be written product over the L 
modes: 



Theorem: consider a (1 + l)-dimensional physical 
theory which depends on a set of real parameters g = 
(g 1 ,g 2 ,...), such that 

• there is a non-trivial conformal point g* , for which 
the model is conformally invariant 

• the deformations from g* in parameter space in the 
positive direction of a given unitary vector e pre- 
serve part of the conformal structure of the model, 
i.e., the eigenvalues of the reduced density matrices 
of the vacuum p(gi) are still of the form given in 
cq.(3.2) for values of the parameters gi ~ g* + ae 



L-l 



* £ =ii(i 



+ e" 



fe=0 



(3.11) 



V^Cg) 



< 0, where q(g) are the correspond- 



yi 



ing parameter-dependent conformal g-factors. 



Once the density matrix of the subsystem is well char- 
acterized with respect to its size L, it is not difficult to 
prove that pl -< pv if L > V . In order to see this, 
we will fix the attention in the majorization within each 
mode and then we will apply the direct product lemma 
from appendix A for the whole subsystem. We initially 
have to observe the behavior in L of the biggest probabil- 
ity defined by each individual distribution for each one 
of the modes, that is, P£ = 1/Z% = (1 + e -6 ") -1 , for 
a = 0,1, . . . , L — 1. It is straightforward to see that 



Then, away from the conformal point there is contin- 
uous majorization of the eigenvalues of the reduced den- 
sity matrices of the ground state along the flow in the 
parameters g in the positive direction of e, i.e., 



9i =9 + ae, g 2 



pG?iHp(<?2), , . 

+ a e, a > a. 



dPl 
dL 



de a 



(1 



dL 



< 



(3.12) 



which implies that P£ decreases if L increases Va. This 
involves majorization within each mode a = 0, 1, . . . , L — 
2 when decreasing L by one unit. In addition, we need to 
see what happens with the last mode a = L — 1 when the 
size of the system is reduced from L to L—l. Because this 
mode disappears for the system of size L—l, its prob- 
ability distribution turns out to be represented by the 
probability vector (1,0), which majorizes any probabil- 
ity distribution of two components. Combining these re- 
sults with the direct product lemma from appendix A, we 
conclude that this example for the quantum XX-modcl 
provides a similar situation for a model with a boundary 
to the one presented in our previous theorem. 



Proof. - 

If the eigenvalues are assumed to be of the form 
given by eq.(3.2), then it is straightforward to see that 



(v ? Ax(5) 



> 0, which assures that the first cu- 



mulant fulfills majorization. The rest of the analysis is 
completely equivalent to the one presented in the previ- 
ous proof of the theorem in sec. Ill, which also proves this 
theorem. □ 

The applicability of this theorem is based on the condi- 
tions we had to assume as hypothesis. Indeed, these con- 
ditions are naturally fulfilled by many interesting models. 
We now wish to illustrate this point with the analytical 
examples of similar situations for the Heisenberg and XY 
quantum spin chains with a boundary. 
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A. Analytical majorization along the anisotropy 
flow for the Heisenberg quantum spin chain with a 
boundary 

Consider the Hamiltonian of the Heisenberg quantum 
spin chain with a boundary 



(4.2) 



where A is the anisotropy parameter. This model is non- 
critical for A > 1 and critical at A = 1. From the pure 
ground state of the system, it is traced out half of it, 
getting an infinite-dimensional density matrix which de- 
scribes half of the system (N/2 contiguous spins in the 
limit TV — > oo). The resulting reduced density matrix 
Pa can be written as a thermal density matrix of free 
fermions [9] , in such a way that its eigenvalues are given 
by 



Pa ("0,^1, • ■ ■ ,«oo) 



1 



with dispersion relation 



2k arcosh(A) , 



(4.3) 



(4.4) 



and rife = 0, 1, for k = 0, 1, . . . , oo. The physical branch 
of the function arcosh(A) is defined for A > 1 and is a 
monotonic increasing function as A increases. On top, 
the whole partition function can be decomposed as 
an infinite direct product of the different free fermionic 
modes. 

From the last equations, it is not difficult to see that 
Pa -< PA' if A < A'. Fixing the attention in a particular 
mode a, we evaluate the derivative of the biggest proba- 
bility for this mode, = (1 + e~ Ca )~ 1 . This derivative 
is seen to be 



dP 



A 



2a 



dA (l 



.) 2 VA2 



1 



> 



(4.5) 



for a = 1, 2, ... oo and for a = 0. It follows from this 
fact that all the modes independently majorize their re- 
spective probability distributions as A increases, with the 
peculiarity that the Oth mode remains unchanged along 
the flow, being its probability distribution always (^, ^). 
The particular behavior of this mode is the responsible 
for the appearance of the "cat" state that is the ground 
state for large values of A (in that limit, the model cor- 
responds to the quantum Ising model without magnetic 
field). These results, together with the direct product 
lemma from appendix A, make this example obey ma- 
jorization along the flow in the parameter. 



spin chain. Let us consider the quantum XK-model with 
a boundary, as described by the Hamiltonian 

oo 

H = - J2 (t 1 + + 0- - 7)«+i + 2A<) , 

n=l 

(4.6) 

where 7 can be regarded as the anisotropy parameter 
and A as the magnetic field. The phase diagram of this 
model is shown in fig.(l), where it is seen that there exist 
different critical regions depending on the values of the 
parameters. Consider the ground state of this Hamilto- 
nian of infinite number of spins, and trace out half of the 
system (if the size of the system is N, we trace out N/2 
contiguous spins, and take the limit N — * 00), for given 
values of A and 7. The resulting density matrix P(\ n ) 
can be written as a thermal state of free fermions, and 
its eigenvalues are given by (see [9]): 



P(A, 7 )("-0,«1, ■ ■ ■ >™oo) = 



1 r.cc 



J (A,7) 



, (4-7) 



where ~ 0, 1, and the single- mode energies e& are 
given by 



<7 



2ke , if A < 1 

{2k + l)e , if A > 1 



(4.8) 



with k = 0, 1, . . . , 00. The parameter e is defined by the 
relation 



e = 7r 



J(VT^a?) 
I(x) 



(4.9) 



I(x) being the complete elliptic integral of the first kind 

I(x) = / (4.10) 
Jo ] Jl-x 2 sm 2 (0) 

and x being given by 

_ /(VA 2 +7 2 -~l)/7 , if A < 1 
X (7/(v / A 2 +7 2 -l) , if A > 1 , 



(4.11) 



with the condition A 2 + 7 2 > 1 (external region of the 
BM-circle [11]). 

We note that the probability distribution defined by 
the eigenvalues of P(\^) is the direct product of distribu- 
tions for each one of the separate modes. Therefore, in 
order to study majorization we can focus separately on 
each one of these modes, in the same way as we already 
did in the previous examples. We wish now to consider 
our analysis in terms of the flows with respect to the 
magnetic field A and with respect to the anisotropy 7 in 
a separate way. 



B. Analytical majorization along uniparametric 
flows for the quantum XY-model with a boundary 



1. Flow along the magnetic field X 



Similar results to the one obtained for the Heisenberg 
model can be obtained as well for a more generic quantum 



We consider in this subsection a fixed value of 7 while 
the value of A changes, always fulfilling the condition 



6 



1 -i 



XX Ising 
critical XY 




critical Ising 



i y 

Figure 1: phase diagram of the quantum XY -model. 



A +7 > 1. Therefore, at this point we can drop 7 from 
our notation. We separate the analysis of majorization 
for the regions 1 < A < 00 and +yl — 7 2 < A < 1 for 
reasons that will become clearer during the example but 
that already can be realized just by looking at the phase 
space structure in fig.(l). 

a. 1 < A < 00.- We show that p\ -< py if A < A'. 
In this region of parameter space, the biggest probability 



b. +yl — 7 2 < A < 1.- For this case, we show that 
Px -< P\' if A > A'. In particular, the probability distri- 
bution for the Oth fermionic mode remains constant and 
equal to (5, 5), which brings a "cat" state for low values 
of A. Similar to the latter case, the biggest probability 
for mode a is P" = (1 + e~ e °) -1 , with 



2an y — , — '- = 2ae 
I(x) 



(4.15) 



and x = {■sj A 2 + j 2 — l)/7. Its derivative with respect 
to A is 



dPl 
dA 



2ae 



de 



(! + e -2ae)2 dA 



(4.16) 



It is easy to see that this time (4r ) < 0, and therefore 

< for a = 1,2,..., 00, which brings majoriza- 
tion individually for each one of these modes when A 

decreases. The mode a = needs of special attention, 

dP?=° 



dP£ 
dX 



from eq.(4.16) it is seen that 



dA 



0, therefore the 



probability distribution for this mode remains constant 
and equal to (|, 5) all along the flow. This is a marginal 
mode that brings the system to a "cat" state that ap- 
pears as ground state of the system for low values of A. 
Notice that this peculiarity is rooted on the particular 
form of the dispersion relation given in equation (4.8), 
which is proportional to 2k instead of 2k + 1 for this re- 
gion in parameter space. These results, together with the 
direct product lemma from appendix A, prove that this 
example fulfills also majorization. 



for the mode a is P" = (1 



) 



with 



(2a + 1)„WLJ?1 = (2a + l) t 
I(x) 



(4.12) 



where x = 7/ ( \/ \ 2 + 7 2 — 1) . The variation of the 
biggest eigenvalue with respect to A is 



dP« (2a + i) e -(2"+i)e de 



dA 



It is easy to see that 



-(2a+l)e^ 2 dA 



(4.13) 



de 
dA 



de da; 
dx dA 

" / dJ(\/T 



x 2 ) 



i(x) y 



dx 



> 




since both (^j§) < and (ff) < 0. Therefore, ^ > 
for a = 0, 1, ... , 00. This derivation shows mode-by- 
mode majorization when A increases. Combining this 
result with the direct product lemma from appendix A, 
we see that this example obeys majorization. 



2. Flow along the anisotropy 7 

In this subsection, the magnetic field A is fixed and 
the anisotropy 7 is the only free parameter of the model, 
always fulfilling A 2 + 7 2 > 1. Thus, at this point we can 
drop A from our notation. We will see that p 7 -< py if 
7^7'i m the two regions 1 < A < 00 and — 7 2 < 

A < 1. In particular, in the region +yl — 7 2 < A < 1, 
the probability distribution for the Oth fermionic mode 
remains constant and equal to (|, J). Let us consider the 



biggest probability for the mode a, P" = (1 
with e a = ue, where 

if A < 1 
if A > 1 , 




(4.17) 



and e as defined in the preceding sections. It is easy to 
verify that 



dP« 



toe 



de da; 



d7 (1 + e-^") 2 dx d7 



< 



(4.18) 



for a = 0, 1, . . . , 00 if A > 1 and for a = 1, 2, . . . , 00 if 
A < 1. The mode a = for A < 1 needs of special atten- 



tion: it is seen that 



dp? 



dX 



0, therefore the probability 



7 



distribution for this mode remains constant and equal to 
(5,5) all along the flow. These results, together with 
the direct product lemma from appendix A, show that 
this case also obeys majorization along the flow in the 
parameter. 



V. CONCLUSIONS 

In this paper we have provided in a rigorous way ex- 
plicit proves for all the majorization conjectures raised 
by Latorre, Liitkcn, Rico, Vidal and Kitaev in previous 
papers on quantum spin chains [5, 6]. In particular, we 
have developed a completely general proof of majoriza- 
tion relations underlying the structure of the vacuum 
with respect to the size of the block L for all possible 
(1 + l)-dimensional conformal field theories. An example 
of a similar situation has been given with the particu- 
lar case of the XX-vaode\ with a boundary, for which the 
explicit calculation of the eigenvalues of the reduced den- 
sity matrix can be performed. We have proven as well the 
existence of a fine-grained entanglement loss for (1 + 1)- 
dimcnsional quantum systems along uniparamctric flows, 
regarded that perturbations in parameter space preserve 
part of the conformal structure of the partition function, 
and some monotonicity conditions hold as well. Again ex- 
amples of similar situations have been provided by means 
of the Heisenberg and XY models with a boundary. Our 
results provide solid mathematical grounds for the ex- 
istence of majorization relations along RG-flows under- 
lying the structure of the vacuum of (l+l)-dimcnsional 
quantum spin chains. 

Understanding the entanglement structure of the vac- 
uum of (1 + l)-dimcnsional models is a major task in 
quantum information science. For instance, spin chains 
like the ones described in the particular examples of this 
paper can be used as possible approximations to the com- 
plicated interactions that take place in the register of a 
quantum computer [12]. Entanglement across a quan- 
tum phase transition has also an important role in quan- 
tum algorithm design, and in particular in quantum al- 
gorithms by adiabatic evolution [13]. On top, the prop- 
erties of quantum state transmission through spin chains 
are also intimately related to the entanglement proper- 
ties present in the chain [14]. Consequently, our precise 
characterization of entanglement in terms of majorization 
relations should be helpful for the design of more power- 
ful quantum algorithms and quantum state transmission 
protocols. 

It would also be of interest trying to relate the results 
presented in this paper to possible extensions of the c- 
theorem [8] to systems with more than (l+l)-dimensions. 
While other approaches are also possible [15], majoriza- 
tion may be a unique tool in order to envisage irreversibil- 
ity of RG-flows in terms of properties of the vacuum 
only, and some numerical results in this direction have 
already been observed in systems of different dimension- 



ality along uniparamctric flows [16]. New strict mathe- 
matical results could probably be achieved in these situa- 
tions following the ideas that we have presented all along 
this work. 
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Note added: after completion of this paper a similar 
work appeared [17] in which entanglement and majoriza- 
tion are considered from the point of view of boundary 
conformal field theory, and where it was noticed that 
there was an error in the proof of the first theorem of 
a previous version of this paper. That error has been 
corrected in this new version. 



Appendix A: LEMMAS ON MAJORIZATION 

This appendix includes the formal definitions of ma- 
jorization [4] as well as a lemma that is used along the 
examples presented in this work. 



1. Definitions 

Let x, y G R N be two vectors such that Yli=i x i = 
Si=i Hi = 1) which represent two different probability 
distributions. We say that distribution y majorizes dis- 
tribution x, written x -< y , if and only if there exist a set 
of permutation matrices {Pj} and probabilities pj > 0, 
Y^Pj = 1, such that 

Y^l'.l'.si- (Al) 

3 

Since, from the previous definition, x can be obtained by 
means of a probabilistic combination of permutations of 
y, we get the intuitive notion that distribution x is more 
disordered than y. 

Notice that in (Al), J2jPjPj = D defines a doubly 
stochastic matrix, i.e. D has nonnegative entries and 
each row and column sums to unity. Then, x -< y if and 
only if x = Dy, D being a doubly stochastic matrix. 

Another equivalent definition of majorization can be 
stated in terms of a set of inequalities between the two 
distributions. Consider the components of the two vec- 
tors sorted in decreasing order, written as (zi, . . . z/v) = 
zK Then, tit* -< y* if and only if 
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k k 

^Xi<Y,yi k = l...N . (A2) 

i=l i=l 

All along this work, these probability sums are called 
cumulants. 

A powerful relation between majorization and any 
convex function / over the set of probability vectors 
states that x -< y => f (x) < f (y). From this re- 
lation it follows that the common Shannon entropy 
H(x) = — ^2i=i x i log 2 ^ °f a probability distribution sat- 
isfies H (x) > H (y) whenever x -< y. In what follows we 
present a lemma that is used all along our work in the 
different examples that we analyze. 



2. Direct product lemma [5] 

If pi -< p 2 , qi -< q_2 then (pi ® gi) -< (p 2 <8> (fa)- This 
means that majorization is preserved under the direct 
product operation. 

Proof. - 

If pi -< f>2 and q\ -< q 2 then p\ = D p f>2 and q\ = 
D q q2 where D p , D q are both doubly stochastic matrices. 
Therefore (pi<£><fi) = {D p ®D q )(p2®q2), where (D p ®D q ) 
is a doubly stochastic matrix in the direct product space, 
and so (pi <X> q\) -< (p 2 ® q 2 )- □ 
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